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Abs t ract 


The He l t 1 e r-London curves for the Interaction of two hydrogen 

* 

atoms are calculated as a function of the Internuclear distance and 
velocity. Electron exchange terms are found to dominate the Inter- 
action energy at laboratory energy 1 kev, to be quite significant 
at 25 kev, and to be negligible at 400 kev. The exchange energy Is 
attractive at 25 kev and below, but at 100 kev It Is repulsive. 
Hence In going from 25 to 100 kev the Singlet and Triplet curves 
go from attractive and repulsive to repulsive and attractive, 
respectively, suggesting that an oscillatory behavior of the cross- 
section versus energy is possible. These results have Important 
Implications for the calculation of atom-atom excitation cross- 
sections In this energy range. 


I nt roduct i on 


Recently there have been a number of experimental and 
2 

theoretical papers on lOn-atom collisions and Lyman- p ro- 
duction and polarization. Activity has also begun in the 

7 k C 

understanding of neu t ra 1 -neu t ra 1 collisions ’ and the theory^ 
of heavy particle scattering in general. Theoretically both 
sysjtems are subject to the difficulty that for a larger than 
two-electron system the electronic wave function at zero ve- 
locity of the nuclei is either unknown (the wave function for 

+ * 

H -H is known exactly and to a degree of accuracy which is 

virtually exact for H-H^) or known only in the interior, 

molecular binding regions for certain selected molecular 

symmetries (for three or more heavy bodies) to degrees of 

accuracy which may permit the calculation of elastic or target 

excitation of rota t i on- v i b ra t i on states scattering but which 

do not permit the study of chemical reactions^, and unknown 

entirely in the region of intermediate Internuclear separation. 

Added to these difficulties of wave function ignorance is the 

fact that adiabatic state* calculated in a molecule fixed frame 

are not correct with respect to the center of mass frame at 

8 * 

infinity and must be transformed into the correct asymtotic 

9 

states . 

An intermediate approach to the heavy particle scattering 
problem would be to divide the molecular space into the interior 
region for which accurate molecular states would be required 
and the intermediate-asymtot 1 c region whose physics can be 
described by use of the states of the separated systems which 


k 


are properly an 1 1 -symmet r 1 zed to account for the Identity of 
the electrons, which becomes Important In the region of small 
overlap, I.e. the unknown "intermediate" region. This is just 
the He 1 1 1 er-London ' ® wave function whose exchange property 
permitted the first qualitative calculation of the binding 
energy of the hydrogen mo.lecule. Once they become available 
the accurate. Interior functions can be fitted to the Heitler- 
London functions in the intermediate region. Since the Heitler- 
London functions are the correct states at infinity such a 

fitting would be equivalent to the transformation of molecular 

~ 9 

states into correct asymtotic states at infinity , and in 

addition since the fitting cOuld be performed before infinity, 

while the interaction was still on but weak, the physics of this 

region would be known. 


Theory and Numerical P rocedures 


Recently Flannery 1 ' has calculated H-H Is1s-1s2s excitation 
cross-sections in the coupled atomic state impact parameter pic- 
ture. Electron exchange is not considered, and the set of direct 
product states used^are precisely those of the van der Waals 
region in which there is no overlap. Qualitatively one would 
expect overlap to be important at close enough encounters to 
produce an excitation of one of the atoms. The justification 
usually given is that at high proton velocities the time needed 
for the electrons to become nonlocal ized in the two-center, 
molecular space is not available. In this paper we calculate 
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the diagonal exchange matrix element, 

j. 

Hi( PtC ls(lA)h/?j)|(R''-^ _ ix" 1 rx," 1 )\h(ie)h(zA) > 



where \/ is the relative velocity vector, of magnitude i the laboratory 
velocity, and protons A and B are going in the directions of negative 
and positive velocities respectively. For mathematical simplicity 
we consider the special case of zero impact parameter (b»0, where 
R*(b +(vt) )* In the rectilinear trajectory picture) cor respond! ng 
to a head-on collision. This avoids the m coupling in the integral 
over the product containing r^ exp(lv/r_), because v is taken along 
the polar axis. Then we can write, 


= '2(;) l -(2k*') A(4f) 

M 


where Jk(vR/2j) is the spherical Bessel function and Pj(*|)the Legendre 
polynomial.The overlap and nuclear attraction integrals are straightforwardly 
Integrated numerically; and a check with the known results^ for 
zero velocity shows agreement to six decimal places. The electron- 
electron repulsion integrals are more difficult. Use of the Neu- 
mann expansion 1 ^ for r~]| results in an infinite series In partial 
1%, "waves" because of the coupling of the Legendre polynomials 



In the Neumann expansion with those In the plane wave expansion 
(2). It was found that for v (laboratory velocity) -.2 (for 
v-0 only the k-0 and k-2 waves are nonvan l sh I ng , and the 
numerical Integration shows agreement with the exact analytic 
result 1 *^ to five places) the Integrals converged to better than 
three places after four waves (k-0, 1, 2, 3)* For v«l convergence 
was obtained to better than three places after 5 waves, and for 

v*»2 after 6 waves. For v» . 2 the k-1 wave falls off by two orders 

¥ » 

of magnitude from the k-0 wave, and the k*2 falls off by three 
orders. For v-1 each wave after k-0 falls off by one order. 
Convergence was much slower for v-2, with waves after k-0 falling 
off by about i . For v-4 the Integral peaked out at k-1 and 
gradually fell off. Except for R-1.0 the v-2 results are good to 
better than two places, and the v-4 results are good to only about 
one place. All other finite velocity results are good to better 
than three places. 

The He l t 1 e r-London Interaction at zero velocity is given 

by, 

11.(0 ± 

i iTcir 7 



whe re , 
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At nonzero velocity it is given by, 

R.lO t 

TTsi^v) y 

V 

Hq(R) has been recalculated by Flannery and Levy'^ and is the 
first diagonal matrix element in the coupled state calculation 
of Flannery'*. In Fig. 1 Hg(R) and Hj(R,v) are plotted as 
functions of R for laboratory velocities, 0, .2, t, 2, and 4. 

In Fig. 2 Hq(R) + Hj(R,v) and Hg(R) - Hj(R,v) singlet and triplet 
i nteract i Oi.s respectively are plotted for the same values of v. 
The numerical values for the individual matrix elements are given 
in Tables 1-5. 


Discussion of Results 


The most important aspect about the curves in Fig. 1 is 
the fact that the velocity dependent phase factors do not 
effectively reduce Hj until quite high laboratory velocities. 
The trend of the curves is what would be expected on physical 
grounds. Greater velocities permit greater penetration of 
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the atoms, reflected in the shifted minima to smaller R. At 100 kev 
(v i 2} the exchange interaction has no negative (attractive) character, 

9 

except above R » 2.5, where it is small and negative and approaches the 
axis from below. At 400 kev (v » 4) exchange is negligible. Fig. 2 
shows the interesting result that at v ■ 1 and below the exchange forces 
show the familiar attractive Singlet well and the repulsive Triplet 
curve, whereas at v ■ 2 the Triplet curve is now attractive and the 
Singlet repulsive, suggesting the possibility of an oscillatory structure 
of the cross-section versus energy, observed in proton-helium direct 
excitation cross-sections 1 ^ to the higher excited 1 P and ” 4 ) states and 
interpretable in the light of the above possibility as competing Singlet- 
Triplet scattering from the virtual charge transfer channel. Also of 
interest is the crossing of the T(v ■ l) and S(v ■ l) curve at about 
R * 2.4. Calculation of the excited state matrix elements for H-H and the 
charge transfer matrix elements for H + -He and the excitation cross-sections 
for these systems will answer this question and will be the subjects of 
later papers. It is clear from this study that electron exchange, long 
understood as the force which binds the hydrogen molecule and known to be 
Important in low energy scattering, is also important in atomic collisions 
in which the nuclei and the bound electrons have comparable velocities. 
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Figure 1 


Figure 2 


Figure Capt l ons 


. Hq(R) and Hj(R,v) versus internuclear distance (atomic 
units). Hj curves are labelled for laboratory velocities, 
0, .2, 1, 2, bnd k . 

. T (Triplet) and S (Singlet) resultant curves versus 

internuclear distance for the laboratory velocities 

Indicated. The small hump in the T^V-O^and T(V = *2^are 
the result of the drawing. 
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